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Abstract. We measure the flow of water through mixed packings of glass spheres and soft swellable hydrogel 
grains, at constant sample volume. Permeability values are obtained at constant sample volume and at 
porosities smaller than random close packing, for different glass bead diameters D and for variable gel 
grain diameter d, as controlled by the salinity of the water. The gel content is also varied. We find that the 
permeability decays exponentially in n(Z)/d)' where n = Ngei/Ngiass is the gel to glass bead number ratio 
and b is approximately 3. Therefore, flow properties are determined by the volume fraction of gel beads. 
A simple model based on the porosity of overlapping spheres is used to account for these observations. 
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1 Introduction 



The flow of water through packings of granular materi- 
als where swellable or soft grains are added is a question 
relevant to a wide range of problems. In the field, the pre- 
diction and modeling of the flows through natural soils 
that contain clay, which swells with water, is still sub- 
ject to active research by geophysicists |1I2| . On the other 
hand, the control of both the drainage and the water re- 
tention capability of granular materials where hydrogels 
have been added has recently been proposed. The appli- 
cations of such additives span a wide range, from baby 
diapers where polyacrylamide gels are added to the flber 
pads [3] , to the sandy soils of desertifled countries where 
hydrogels have been envisioned as a way to reduce the 
volume of irrigation (for a review, see |4I5) ). In the latter 
case, the efficacy of superabsorbents as soil additives was 
shown to strongly depend from one trial to another and 
is still subject to controversy. 

The properties of flows through packings of hard grains 
have been extensively studied, both theoretically [6 7 8', 
numerically [9110111] and experimentally [12113] . Based on 
the pioneering work of Darcy on the bulk resistance of a 
porous medium to incompressible viscous fluid flow, the 
superficial fluid velocity U = Q/S, where Q is the vol- 
ume flow rate across a section S of granular medium, is 
related to the pressure gradient VP = AP/Al in the flow 
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direction by the following linear equation: 
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KD^ 



(1) 



Here rj is the dynamic viscosity of the fiuid, k is the per- 
meability, with the dimension of an area, and K is the 
intrinsic permeability, a dimensionless number character- 
izing the packing properties. The term accounts for 
the dependence with the inversed squared specific surface 
of solids in Darcy's law, and D is typically a particle size. 
Hence, for the simplest case of spherical grains, D is taken 
as the mean diameter [14| . 

Much research has been devoted to understanding the 
value of the permeability. For a random close packing of 
spheres, where the porosity is e ~ 0.36, the comprehensive 
experimental study reported in Ref.114] gives: 



Ko = 6.3 X 10" 



(2) 



From here, the dependence of the permeability on the 
porosity e, the average pore size, and the tortuosity was in- 
vestigated for grains of various shapes: spheres [1511419116111] , 
aspherical particles [T^, sand [T3] or rocks and sintered 
glass beads [T7]- Several quantitative relations were pro- 
posed and tested, mostly for porosities larger than the ran- 
dom close packing case. A review of these models can be 
found in [18) . Among them, the most known and used was 
established by Kozeny [6] and Carman fTl and relates the 
intrinsic permeability K to the porosity e. The Kozeny- 
Carman model proposes the following relationship: 
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In this model, the pore space is seen as a bundle of paths 
which do not intersect each other. This is of course unreal- 
istic since in most porous media, the pore space is highly 
interconnected. Still, it is in fairly good agreement with 
experimental studies for porosities close to random close 
packing, as in Refs. |12in] and, in particular, for e = 0.36, 
Kkc in Eq. [3] agrees well with Kq in Eq. [21 More recently, 
the percolation theory [lO^ and self-similar modelling j20i 
were applied to model the porous media as an attempt to 
get more realistic descriptions. 

When soft and/or swellable grains are added the key 
question of the porous volume becomes challenging. In- 
deed, the shape and volume of soft grains adjust to the 
load, making possible a porosity value smaller than the 
random close packing one, down to the full clogging case 
[H] . Relevant to us are the works of Weissberg [T5] , and 
of Martys et al. [5] where packings of overlapping spheres 
were considered, using statistical calculations in the for- 
mer case, and simulation in the latter one, yielding porosi- 
ties as small as 0.1. An important result of [15] is that the 
permeability of overlapping spheres of increasing volume 
is expected to decrease exponentially with the volume of 
the spheres. Also of interest to us is the case of sintered 
glass beads studied experimentally and theoretically by 
Wong et al. [T^- In this work, pores are seen as tubes of 
random size that are interconnected, and the reduction of 
the porous space is obtained by systematically reducing 
the radius of a random tube element by a fixed factor. 
The result is a power-law dependence of the permeability 
with porosity, but with an exponent that increases towards 
smaller porosities in an unknown way. The main result 
however is that the scaling behavior of the permeability 
of sintered glass bead samples and rocks is determined by 
the skewness of the pore-size distribution. The similarity 
between the two types of samples arises from the forma- 
tion process, which reduces more the larger pores that the 
smaller ones. 

It is unclear whether or not either of the above mod- 
els, overlapping spheres or systematically reduced pores, 
is applicable to our system of mixtures of hard and soft 
swellable grains. Unfortunately, very few experiments have 
been performed on such systems, and no systematic study 
can be found in the literature where every parameter was 
carefully controlled. Measurements of the permeability of 
sandy soils containing hydrogels were reported by Bhard- 
waj et al. [5]. The authors qualitatively show that the 
swelling, the size and the amount of gel added play a 
role on the measured permeability. However, the swelling, 
tuned by the ionic content of the draining solution, wasn't 
controlled during the course of the experiment although 
hydrogels release ions at the beginning of the drainage. In 
addition, mixtures of gels and sand tend to expand upon 
swelling of the gels with water. The variations in volume of 
the sample appear to be strongly history dependent when 
the mixture is in a cylindrical column, due to wall effects. 
This varied volume impairs the analysis because both the 
porous volume and the solid (gel) volume increase. 
In this paper, to circumvent the above difficulties, we de- 



scribe permeability measurements with a controlled sys- 
tematic loading and swelling procedure. In particular, novel 
features include: (a) The volume of the pile is forced to 
remain constant, dry or wet; (b) First column filling is per- 
formed under vacuum to ensure saturation of the porous 
volume, and a highly salty water is first used; (c) Then, 
a gradual swelling of the gel is induced by flushing in 
less salty water; (d) Ionic contents of both the feeding 
solution and the collected solution after the pile are con- 
stantly monitored; (e) Rinsing with a new feeding solution 
is performed until the ionic content at the outlet reaches 
the feeding value; (f) Degased salt solutions are used, and 
stored in a reservoir where nitrogen is bubbled to avoid 
CO2 bubbles nucleation. By systematically varying the 
amount of gel added, their swelling and the size ratio of 
soft and hard grains, we gradually modified the porosity 
of our mixtures, and were able to determine a functional 
form to describe the decrease in permeability upon soft 
(gel) particles addition from the value obtained with the 
initial packing of spherical hard (glass) beads. 



2 Experimental set-up 

We use a model system of spherical glass beads and com- 
mercial hydrogels. The mixture is saturated with water 
of controlled salinity to set the osmotic pressure that, to- 
gether with the mechanical pressure, controls the swelling 
of the gels [22123] . The gels are forced to swell within the 
porous space between the densely packed hard spheres: 
this way, we probe the clogging produced by the soft par- 
ticles. The experiment (Fig. [T]) consists in pushing wa- 
ter (with salt) through a water-saturated mixture of glass 
beads and gels contained in a cylinder at a controlled flow 
rate, while measuring the pressure drop over a given height 
difference along the column length. The experimental de- 
tails are discussed below. 



2.1 Materials 

The glass beads (Potters Industries, PA) are spherical. 
Three bead diameters D were used: 0.5, 1, and 2 mm. 
The size distributions are measured using an automated 
sizer (CamSizer) and fit to a Gaussian curve with a width 
of order 100 /im and a mean diameter D. We also mea- 
sured the average mass of one glass bead for each size 
by weighing several hundreds of carefully counted beads. 
Prior to being used, the glass beads are treated so as to 
fully hydrate their surface. They are first burnt in a fur- 
nace at 500C for 72 hours and then acid-washed in a IM 
HCl bath for 1 hour. The beads are then thoroughly rinsed 
with pure deionized water and dried in a vacuum oven at 
HOC. 

The gel particles were kindly provided by Degussa Inc. 
(Stockosorb SW) and are shown in Figure [TJ They con- 
sist of a cross-linked copolymer of potassium acrylate and 
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Fig. 1. (Left) Photograph of the experimental system: 
dry gels (ddri/=0.2 mm), swollen gels (d=2 mm) and glass 
beads (D=0.5 mm). (Right) Schematic representation of 
the experimental set-up. A cylindrical column holds the 
mixture of gels and glass beads. Volume is kept constant 
by filling the extra space on top with a sponge (1-cm high) 
and centimer-sized balls. Water is pumped into the top of 
the column at a constant volume flow rate Q. Pressure 
transducer measures the pressure drop AP along a Al = 
50.8 mm height section of the sample. 
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Fig. 2. (a) Volume of the swollen hydrogels as a function 
of KCI concentration in water. Right axis: swelling ra- 
tio. At lower salt content, the swelling ratio plateaus. At 
higher salt content, it fits to a power law (line) with an ex- 
ponent -0.46. Dotted line: theoretical prediciton from [22 . 
(b) Mixtures of 0.5 mm glass beads and hydrogels at var- 
ious gel to bead number ratios: Normalized permeability 
K/Kq as a function of KCI concentration in the draining 
water. The dashed lines indicate the plateau value at low 
salt content. 



acrylamide (50/50). As a result of the manufacturing pro- 
cess that produces the particles by grinding a bulk gel and 
then sieving, their shape is faceted and random. The dry 
particles are further sieved between 212 and 300-micrometer 
sized meshes. The swelling ratio of these particles is mea- 
sured for KCI solutions of increasing salt concentrations 
as follows. Batches of 100 particles were carefully counted, 
then soaked into a solution of given salinity in a tared 
weighing dish. After an hour, the excess of liquid was 
sucked up with a pipette and a paper wipe. The mass 
of the swollen particles is then measured with a precision 
scale. The process is repeated for every (higher) salt con- 
centrations. Three different batches of 100 particles were 
used and we eventually obtain the average mass of one gel 
particle at increasing salt contents. We also measured the 
average mass of one dry gel particle for over 300 speci- 
mens: m — 14.8 ± 0.5 fig. The swelling ratio of the gel is 
defined as the ratio of the swollen to the dry mass. The 
average volume Vgei is calculated from the mass using the 
density of the salt solution as a rough approximate for the 
particle density. The results are shown in Figure [2]-a) . At 
low salt concentration, the swelling ratio is insensitive to 
the salt concentration and plateaus: the cross-linked poly- 
mer network is swollen to its maximum extent. For [KCL] 
values larger than 2 x 10~^ mol.L^^, it decreases and the 
data can be fit to a power law with an exponent -0.46. 
This result is in good agreement with other experimental 
data reported in the literature [35] for high salt concentra- 
tions that yield an exponent ranging between -0.3 (dotted 
line, theoretical value) and -0.6. It can be noted that in 
the range of salt content used here the volume of the gel 
particles changes by one order of magnitude. 



The permeability measurement apparatus, shown in Fig- 
ure [TJ consists in a custom-made cylindrical column in 
polycarbonate (height 30.5 cm, inner diameter 7 cm, McMaster- 
Carr). It is closed at the bottom end by a brass sieve 
(mesh size 90 micrometers) that is attached to a funnel. 
The samples poured in the experimental column are pre- 
pared by mixing small amounts of glass beads (typically 
25 to 100 g) and dry gel particles to ensure a good homo- 
geneity of the final mixture, to a total bead mass of 700 g. 
Using this bead mass and measuring the volume occupied 
by the dry glass beads in the column, we checked that the 
beads arrange themselves in a statistical packing of poros- 
ity e=0.38±0.04, slightly higher than the usual figure for 
the random close pack porosity of 0.36. Note that this dry 
porosity based on the glass beads, only, is unchanged by 
the addition of even the larger amounts of dry gel. Hence, 
the initial packing of the spheres is not affected by the 
presence of dry gels. For data analysis, we use the average 
masses of individual gels and of individual glass beads to 
convert the gel to glass bead mass fraction that we use in 
practice into a gel to glass bead number ratio n defined 
as: 

n = N gel /N glass (4) 



2.2 Procedures 
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The upper end of the column is sealed by a plug connected 
either to a vacuum pump (filling step) or to a gear pump 
(measurement step). First, the column is filled with a KCl 
solution at high concentration (0.07 mol.L"^) under vac- 
uum and from the bottom. This procedure allows for the 
saturation of the pores with the solution. The solutions 
used are prepared with deionized water, and then degased 
under vaccum to minimize the nucleation of CO2 pock- 
ets within the sample. Nitrogen is bubbled in the water 
reservoir during the entire experiment (filling and mea- 
surement) to prevent the redissolution of CO2 in water. 
During the filling, the volume of the glass bead and gel 
mixture is kept constant by adding a sponge and marbles 
on the top of the sample, and by using a salty solution 
to minimize the gel volume expansion. Consequently, the 
gels swell only within the porous space available between 
the glass beads. 

The upper plug is then connected to a gear pump (Mi- 
cropump Inc.). We use two different pump heads depend- 
ing on the range of flow rates desired: 0.092mL/revolution 
for flow rates Q ranging from 15 to 330 mL/min and 0.017 
mL/min for Q between 3 and 60 mL/min. On the side of 
the column, along its height, two ports are drilled through 
the column wall at a distance Al — 50.8 mm, to allow the 
connection of a wet-wet differential pressure transducer 
(26PCA, Omega). The pressure sensor is calibrated by 
connecting a water column of controlled height. It has a 
typical sensivity of 3.5 mV/kPa. The range of Q is chosen 
so as to probe a pressure drop variation of one order of 
magnitude. We carefully checked that the pressure drop 
AP linearly increases with the flow rate Q for every mea- 
surement performed: this shows that all flows are in the 
viscous regime. The pressure drop AP over Al is used to 
compute the permeability k and the intrinsic permeability 
K defined by equation [1] Here, the intrinsic permeability 
K is computed using the glass bead diameter D as the 
typical size that normalizes n. For every change in salin- 
ity of the draining solution, the mixture is rinsed until the 
ionic content at the lower end of the column, measured 
with an electric conductivity-meter, reaches the reservoir 
value. 



3 Results 

First, we focus on the pure glass bead systems as a way to 
test our set-up. The permeability of the glass beads with- 
out gel, Kq, is plotted in Figure[3]for the three D values at 
the abscissa, gel to bead number ratio n=Q. The horizon- 
tal dotted lines indicate the value Kq = 6.3 x 10~^ given 
by Eq|2] Averaged together, our three measurements at 
n = give Kq = (6.1 ± 0.7) x lO"'*, in good agreement 
with Eq. [21 This confirms the integrity of our procedures 
and indicates that the size of the sample is large enough 
to neglect wall and edge effects. 

We now turn to the gel and glass bead mixed systems. 
In a first series of experiments, we used the salinity of the 



draining solution to control the flow properties of a pack- 
ing of 0.5-millimeter glass beads with gels. The series is 
performed with solutions of decreasing salinity. The per- 
meabilities obtained are shown in Figure HJ-b as a function 
of the KCl concentration in the solution. For every gel con- 
tent tested, we observe that the permeability plateaus at 
small salt concentration and then increases. The cross-over 
between these two regimes occurs at a salinity comparable 
to the salinity where the size of a single gel particle starts 
decreasing (Fig.[2]-a). For [KCl] concentrations larger than 
5 X 10"** mol.L^^, the permeability follows the same trend 
but inverted as the swelling ratio of one gel particle. As 
an example, a sample with a mass fraction of 0.4% (or 
number ratio n — 0.06) has a permeability that spans the 
range 2% to 70% when increasing the KCl concentration 
from 10~^ to 1 mol.L^^. This shows that salinity stands as 
an efficient external trigger for the permeability: the typ- 
ical size of the gel particles appears to be a characteristic 
size of the problem. On the other hand, it is important 
to rinse the granular pile before acquiring any data be- 
cause of excess ions in the dry gels. Omitting to take into 
account the excess ions in the dry particles leads to mis- 
interpretations of measurements of the permeability with 
these gels, and to time-dependent measurements as in [5]. 

In addition, we verify the obvious fact that larger gel 
contents yield smaller conductivities. The five curves pre- 
sented in Figure[2]-b correspond to mass fractions of 0.4%, 
0.2%, 0.1%, 0.05% and 0.01% from bottom to top. Let us 
now focus on the low salinity case. This series at constant 
salt is meant to isolate the effect of the size ratio from 
other effects such as the elastic properties of the gels that 
are modified together with the gel size when salt is added 
|22| . The plateau value at small salinity is measured for 
every gel content used and is plotted as a function of the 
gel to glass bead number ratio n in Figure [3]-a. The same 
data sets are collected for two more glass bead sizes (1 and 
2-mm) and are shown in Figs. [2]-b,c. Interestingly, for all 
bead sizes, the permeability decreases exponentially with 
the gel to glass bead number ratio. The data sets fit to an 
exponential decay function with two fitting parameters: 
a characteristic number ratio ric and the permeability at 
zero gel content K^: 

K = Koe-^ (5) 
with Ko = (6.1 ±0.7) X 10--^. 

We use equation [S] to analyze our measurements of per- 
meability as a full function of salinity, for all gel mixtures 
with D — 0.5 mm glass beads. For every salt content, the 
conductivity, normalized by the reference without gel, Kq, 
is plotted as function of n in Figure S) All the data can be 
fit to the exponential decay function given by Eq|S]with a 
characteristic number ratio ric- 

We now turn to the analysis of the characteristic num- 
ber ratio Uc measured in Figures |3] and 21 It was earlier 
found that the gel volume is a characteristic parameter 
of the problem (see Figure [5]) . The two series of experi- 
ments, with different glass bead sizes and with gel sizes 
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Fig. 3. Low salt plateau value of the intrinsic permeability as a function of the gel to glass beads number ratio n 
for three different glass bead diameters: a) 0.5 mm b) 1mm c) 2 mm. All the data can be fit to an exponential decay 
function (full lines) with a characteristic number ratio Uc shown by dashed lines. The dotted line indicates the value 
Ko in Eq. H 
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Fig. 4. Mixtures of 0.5 mm glass beads and gels. Perme- 
ability data from Fig. [5] versus gel to glass number ratio. 
Each curve corresponds to a given salt concentration. Data 
fit to an exponential decay function (lines) with a charac- 
teristic number ratio Uc reported in FigE] as diamonds for 
every KCI concentration. 
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decreased by salt addition, can be seen in terms of a vari- 
ation of the size ratio of the two types of particles. For 
the glass bead size, we take their diameter D while we 
define a typical gel size d at a given salt concentration as 
Vgei = f with Vgei, the swollen volume of gel particles 
taken from Fig.[2]-a. Hence, we plot ric as a function of the 
glass to gel size ratio D/d. We find that the dependency 
of ric with D/d is consistent in both situations, namely 
when this ratio is varied either by decreasing the gel size 
(salt added, diamonds) or by increasing the glass bead size 
(square dots). Moreover, the data can be fit to a power law: 
Tic = a{D/df with a - 0.40 ± 0.04 and 6 - 2.6 ± 0.4 (fuU 
line in Fig. [S|). Interestingly, the exponent h is close to 3. 
In Figure [5l the dotted line corresponds to a power law 
with an exponent 3 and an adjustable prefactor. It nicely 
fits the data within our experimental incertainty with a 
prefactor a* = 0.58 ±0.04. Hence the (D/d)'' dependence 
can be seen as the ratio between the volumes of both types 
of particles, leading to: 



D 



(6) 



Therefore, the argument in the exponential function can 
be recast as a function of the gel to total solid volume 
fraction cb: 



id^ 



glass to gel size ratio D/d 



id^ + 



(7) 



Fig. 5. Characteristic number ratio function of the 

glass bead to gel size ratio D/d. Squares: D is changed, 
d fixed. Tic fits to a{D/d)^ (full line) with an exponent b 
close to 3 and a of order 1. Diamonds: D=0.5 mm, d is 
decreased by KCI addition. All the data fit to a*{D/d)^ 
(dotted line). 



Altogether, we find that the normalized permeability 
K/Kq plotted as a function of the normalized number 
ratio n/ric collapses on a single master curve shown in 
Figure El-a. This master curve can be used to predict the 
drainage reduction provided that the characteristic num- 
ber ratio Tic is calibrated as a function of the glass to gel 
size ratio D/d. The calibration shown on Fig. [5] allows for 
the creation of another master curve where 4> is the unique 
parameter, plotted in Figure [S]-b. 
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gel volume fraction i|)=nd^/(nd^+D^) 



Fig. 6. a) Permeability data for different glass bead sizes and different gel particle sizes collapse on a master curve 
providing that K is normalized by its value witliout gel Kq, and tlie number ratio n by its characteristic value n,. 
obtained from the exponential fit. Full line: exponential decay function (no adjustable parameters). (b) Normalized 

permeability as a function of gel to total solid volume fraction = "( ^^^s'^^ja ) for the three glass bead diameters 
D studied and upon gel shrinking by KCl addition into 0.5-mm sized glass beads systems. Line: exponential decay 

function K / Kq = e^^'fi-*) with a* = 0.58. 



4 Discussion 

Our results show that the permeability of mixtures of soft 
and hard grains can be described by the general law 

K = Kqc -'(o/dY-i ^ (8) 

Here a* — 0.58 ± 0.04 is a numerical constant acounting 
for the gel softness (or equivalently, the salt content), Kq 
is the intrinsic permeability of the glass beads pile given 
by Eq. [2j and (p is the gel to total solid volume fraction. 

In the following, we derive simple expressions to ac- 
count for our results, as summarized by EqlS) Because 
our gels are so soft, we first assume that their effect on 
the packing is to close some pores off while pores that 
aren't affected by a gel particle remain with an unchanged 
conducting capability. Hence, the number of open pores is 
decreased but the size and shape of these open pores re- 
main unchanged. To estimate the new porosity, we use 
the derivation by Weissberg [TS] for overlapping spheres. 
Here we reproduce the main ingredients. The porous vol- 
ume is made of the points in the space that are not con- 
tained by any sphere. The porosity is thus regarded as the 
probability that a random point is not contained by any 
sphere, or, equivalently, as the probability that the por- 
tion of the space which lies within a distance equal to the 
particle radius from the random point contains no sphere 
centers. First, for overlapping, monodispersed spheres of 
volume V, we derive the probability of placing N centers 
at random in a large finite volume V in such a way that 
a smaller volume equal to v contains no centers. In this 
case, each random placement of a center is independent 
of the positions of the other centers, this probability is 
[{V — f)/y]^ = [1 — v/V]'^ . For a larger and larger sys- 
tem where N tends to infinity, holding the sphere volume 
V and the number density N/V fixed, the probability that 
the volume v contains no centers becomes: exp[—Ny]. 



Now, we turn to the two-size problem. We have glass 
beads of size D that are non-overlapping, and gel over- 
lapping particles of volume v =gei= ■^dP ■ The probabil- 
ity for a point to be in the porous space is the product 
of the two probabilities given by each type of particles. 
Because the glass beads don't overlap, the corresponding 
term is simply the initial porosity e. The gel beads are 
assumed to be spherical. The corresponding term is then: 
exphfd^TVge//^]. 

Altogether, we find that the porosity of a composite sys- 
tem with hard glass spheres, and soft, overlapping gel par- 
ticles is given by: e eyi\)[—^d?' Ngei/V]. 
From this, the derivation of the conductivity is made fol- 
lowing the method developped by Carman [7] . Here, we as- 
sume that the pores already present in the glass sphere pile 
are either closed by the gel or left open and unchanged. 
Thus, the conductivity is simply the product of the ini- 
tial conductivity scaled down by the porosity factor, 
yielding: 

K ^ Kqc-^'o"^^^'/^ (9) 

This expression can be recast using the experimental vari- 
ables used to fit our data in FiglHl The total sample vol- 
ume V can be calculated using the glass beads volume 
since dry gels occupy a small enough volume for it to be 
neglected and our experimental procedure prevents the 
sample from swelling. Then, V = NgiassVgiass/i^ - e), 
where Vgiass = is the volume of one glass sphere. 

Hence: 

— CI 1 "gal" gal 
K^Koe ^ ""-glass^gla.. (IQ) 

Our calculation, however rough it is, yields an expo- 
nential decaying function of the total gel volume to the to- 
tal glass volume, in remarkably good agreement with our 
experimental results. The derivation first assumes a spher- 
ical and uniform shape for the gel beads which is quite a 
strong assumption. Second, modelling the gels as overlap- 
ping spheres also assumes a loss of volume when gel is 
squeezed between the glass spheres. Refering to studies of 
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highly swollen gels under compression [23], we know that 
under large strains, gels release solvent and deswell. Here, 
the size of the pores left between glass spheres are about 
ten times smaller than the gel size when unstrained, so the 
strain applied to the gels is high and deswelling probably 
occurs. Also, with our procedure that prevents the pile 
from swelling, the pile is submitted to a load that can be 
as high as needed for the distance between glass beads to 
remain constant. Hence, the effect of an increased osmotic 
pressure upon salt addition is equivalent to an increased 
mechanical pressure yielding deswelling, although salt ad- 
dition leads to both a decrease in size and a increase in 
elastic modulus |23j . This explains why experiments made 
with salty solutions could be analyzed in terms of size 
(FiglS)) without taking into acount the increased stiffness 
of the gel. Moreover, our model assumes that gel particles 
totally close some pores while others remain unchanged. 
In the framework developped by Carman, pores are seen 
as a bundle of parallel tubes of uniform diameter. The 
above assumption means that gels do not partially invade 
a tube. This hypothesis is supported by the large size ra- 
tio between gels and pores and the high softness of the gel 
material (Young's modulus of order 1 to 10 kPa). Hence, a 
gel will swell into any neighboring pore "tube" so as to oc- 
clude it and to get in contact with the next grain. Finally, 
the model neglects excluded volume in assuming that the 
addition of a new gel does not depend on the gel beads 
that are already present. This would increase the slope of 
the decaying exponential behavior. Looking at the argu- 
ment of the exponential function, our experimental results 
yield a prefactor 1/a* = 1.72 while our model predicts a 
prefactor (1 — e) = 0.64. Among other effects, the discrep- 
ancy could in part come from the excluded volume effect 
mentionned above. 



5 Conclusion 

To determine the properties of flows through a mixed 
packing of soft and hard particles, we presented a series 
of measurements of the hydraulic conductivity of a pile 
of glass spheres in which a controlled amount of hydrogel 
particles was added. The volume of the pile was forced to 
remain at the size of the glass bead pile without gels. Care 
was taken to fully saturate the porous volume with water 
at all time, and to precisely control the ionic content of 
the water flowing in the pile. Hence, the efficient volume 
available for conduction and the solid (gel) volume were 
known at all time, thus improving upon prior work. We 
find that the conductivity of the pile is equal to its value 
without gels scaled down by an exponentially decaying 
function of the volume fraction of gel presented in Eq. [S] 
This result can be explained by modelling the gel parti- 
cles as overlapping spheres. The porosity of the assembly 
of hard non-overlapping glass spheres and soft overlapping 
gel particles was derived, and from there, the conductiv- 
ity of the pile was deduced assuming that pores are either 
open or fully occluded by gel. This simple derivation holds 
because our gels are soft and large enough compare to the 
glass beads. 



Turning to the efficacy of gels as drainage modifiers, our 
results allows for the determination of the right amount of 
gel to add to get a given reduction in conductivity. This 
work could be continued by taking into acount the swelling 
of the pile. The experimental device should be carefully 
designed so as to carefully control the volume expansion 
of the pile upon gel swelling, and to make it independent 
of the history of the sample and of wall effects. 
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